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ABSTRACT 

We study the mass distributions, circular velocity curves (CVCs), line-of-sight kine- 
matics and angular momenta out to many R e for a sample of 42 cosmological zoom 
simulations of galaxies with stellar masses from 2.7 x 10 10 AfQ to 4.8 x 10 11 M Q . In 
order to reduce the particle noise at large radii, we temporally smooth the observables 
of the simulated galaxies in a static potential. The mass of the simulated galaxies is 
parametrised by the circular velocity at r = 5 R where R e is the projected stellar half 
mass radius within 10% r va . We find: (i) The projected stellar density distributions at 
large radii can be well fitted by Sersic functions. The Sersic indices range from ~ 3 to 
~ 13 and correlate with stellar mass and galaxy size (low n, low mass, small size), (ii) 
The dark matter halo density profiles are consistent with simple power-law models, 
corresponding to flat dark matter CVCs for lower-mass systems, and rising CVCs for 
high-mass halos. (iii) The massive systems have nearly flat total (luminous plus dark 
matter) CVCs at large radii, while the less massive systems have mildly decreasing 
CVCs. The slope of the circular velocity at large radii correlates with Sersic index and 
circular velocity itself, (iv) The short axes of simulated galaxies and their host dark 
matter are well aligned and their short-to-long axis ratios are correlated, (v) Stellar 
root mean square velocity v Iias (R) profiles are slightly falling, consistent with plane- 
tary nebulae observations in the outer halos of most early- type galaxies (ETGs). There 
are no analogues in the simulated galaxies of the second group of ETGs with rapidly 
falling v rnls (R) profiles, (vi) The line-of-sight velocity fields v show that rotation prop- 
erties at small and large radii are correlated. Most radial profiles for the cumulative 
specific angular momentum parameter A(i?) are nearly flat or slightly rising, with 
values in [0.06, 0.75) from 2R C to 5R C . There are a few cases for which the local 
\v\/a(R) decreases at large R, and some with localalized maxima of specific angular 
momentum. These properties are consistent with observations of ETGs at large radii, 
(vii) Stellar mass, ellipticity at large radii e(5 R G ), and A(5 R e ) are correlated: the more 
massive systems have less angular momentum and are rounder, as for observed ETGs. 

Key words: galaxies: kinematics and dynamics - methods: iV-body simulations - 
methods: numerical 



1 INTRODUCTION 

Recent observations and cosmological simulations suggests a 
two-phase formation scenario for massive early-type galax- 
ies (ETGs), in which an early phase of rapid star forma- 
tion driven by cold accretion and dissipative mergers is fol- 
lowed by a prolonged phase of mass accretion through gas- 
poor major and minor mergers. Observations have shown 
that a population of old, massive (~ 10 n M Q ) and red 



ETGs were already in place at redshifts z = 2 — 3 (e.g . , 
iFontana et"alll200d : HTbert et alj|201(il ; ICassata et al J 120111 ). 
and that these galaxies have sizes several times smaller and 
densities an order of mag nitude higher than present-day 
ETG s of similar mass (e.g - .TPaddi et al.ll2005l ; iTruiillo et all 
120071 ; Ivan Dokkum et all I2OI0I ). Recent simulations have 
found that massive galaxies grow initially through rapid star 
formation fuelled by infall of cold gas at z > 2, leading to 
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an old population of 'in situ' stars. Subsequently, the sim- 
ulated galaxies grow through minor mergers, accreting old 
stars formed in subunits outside the main galaxy halo. The 
accreted stars are preferentially added to the outer haloes 
of the host systems, leading to effic i ent size evolution (e.g. , 
Naab et all 120091; lOser et all |2010|; [Feldmann et al.l l20ld ; 



Hopkins et alj 2010l ~ Johansson et alj 20121 ) . 

While there is relatively little stellar mass in the outer 
regions of ETGs at galactocentric radii > 2 R c , these stars 
may carry indispensable information about the late assem- 
bly history of the galaxies. With relaxation times in the 
present-day outer stellar halos of up to several Gyrs, the 
record of the rec ent halo accretion h i story may be rel atively 
well preserved (|van Dokkuml 120051 ; iDuc et alj l201ll ). The 
outer halo stars may also contain a s ignificant fraction of 
the a ngular momentum of the galaxies (iRomanowskv fc Falll 
2012), and they can be used as gravitational tracers to study 
the mass distribution, dark matter fraction and potential of 
ETGs at large radii. Therefore, a detailed investigation of the 
structure and kinematics of the outer stellar haloes, com- 
paring simulated and observed halos, may reveal important 
information about the formation history of ETGs. 

Traditional long slit measurements and more recent ob- 
servations with integral-field units (IFUs) have provided 
detailed kinematic and dynamical informati on about the 
central regions of ETGs, < 1 — 2R e (e.g.. Bender et al.l 
11994 iGerhard et ai1l200ll ; ICappellari et all 120061 ). Using 2- 



dimensional stellar kinematics from the SAURON IFU out 
to ~ lRc, and a measure of the projected angular mo- 
mentum condensed into the cumulative X(R) parameter, 
ETGs can be separated into two main groups, fast rota- 
tors (FR, MR) > 0-1) and slow rotators (SR, X(R) < 0.1 
lEmsellem et al. I l2007l . l201ll ). This division is part of a wider 
dichotomy between oblate-spheroidal disky, coreless, rotat- 
ing ETGs with little radio and X-ray emission, and triax- 
ial boxy, co red, non-rotating, radio-loud and X-ray bright 
systems (see Bender et alJll989l : iKormendv fe Benderll 19961 ; 
iKormendv et alj|2009h ~ 



In the outer halos, obtaining kinematic information is 
much harder because of the rapid decline of stellar surface 
brightness with radius. Most of the known kinematic prop- 
erties come from observations of Planetary Nebulae (PNe) 
which have been found to be good tra cers of the stars and 
can be observed up to ~ 8R e (e.g., Mendez et alj l200ll ; 
ICoccato etafl 12009 ; iMcNeil-Movlan et al.ll2012j ). Recently, 
individual IFU p ointings and slitlet masks have also been 
used at 3 - 4 R P dWeiimans et al.ll2009l ; iProctor et alll2009l ; 
iMurphv et ai1l201ll ). The PNe observations show that the 
division between FR and SR is largely preserved to large 
radii, although there are also galaxies whose A(_R)-profiles 
drop significantly outwards, possibly implying that the high 
inner values are due to disks whose light contribution fades 
towards large radii. Most of the PNe RMS velocity profiles 
decrease slowly with radius, but a subset of galaxies show 
steep 'quasi-Kepleria n' outer decreases of velocity dispersion 
(|Coccato et alj|2009\ 

Analysis of outer stellar kinematics, strong lensing, and 
hydrostatic equilibrium of X-ray emitting hot gas shows that 
massive elliptical galaxies have nearly isothermal inner mass 



Gerhard et al.l 2001 ; 


Koopmans et al.l 20061; Auger et alJ 


2010l; Churazov et al. 


2008. 20ld; Nagino & Matsushita 



2009). For the lower mass ellipticals, the situation is less 
clear, as the mass-anisotropy degeneracy is stronger for 
declining velocity dispersion profiles, their X-ray emission 
is too faint, and the lensing samples are dominated by 
massive systems. The dynamical modelling of integrated 
light and PNe indicates som ewhat more diffuse dar k mat- 
ter halos in these gal axies dde Lorenzi et alj 120081 . 120091 ; 
iNapolitano et af]|2009l ). Another useful tracer of the outer 
mass distributions is the globular clusters (GCs); espe- 
cially massive ellipticals contain large GC p opulations 
(e.g.. ISchuberth et al.l [20T0I ; IStrader et alj|201ll ). However, 
a larger fraction of GCs may be recently accreted systems, 
as their relative frequency is tilte d more towards sma ll sys- 
tems than that of PNe and light |Coccato et al.|[2012j ). 

The purpose of this paper is to analyze the mass distri- 
butions, outer kinematics, and angular momentum distribu- 
tions of a representative sample of simulated massive galax- 
ies, as a basis for comparison with kinematic observations 
and inferred mass distributions of nearby ETGs. We consider 
a sample of 42 resimulated galaxies from the h igh-resolution 
cosmological simulation of Oser et al] (|2010l ). which grew 
through the two-phase processes of early in-sit u formation 
follow ed by late accretion and minor mergers l|Oser et alj 
2012). The paper is organised as follows: In [J2]we briefly de- 
scribe the present-day (z=0) model galaxies extracted from 
the cosmological zoom simulations, and the method we use 
to derive smooth kinematic maps from their particle distri- 
butions. In [|3] we study the mass density distributions of 
the stars and dark matter in these galaxies, as well as the 
corresponding circular velocity curves (CVCs). Then in Sj?] 
we investigate the observable kinematics in these systems 
out to large radii. We finally consider the cumulative and 
local angular momentum profiles X(R) and \v\/a(R) for the 
stellar components in Sj5] We end by summarizing our results 
in SJ6] 



2 THE SIMULATED GALAXIES AND THEIR 
KINEMATIC ANALYSIS 

The "galaxies" studied here are extrac t ed from the cosmo- 
logical zoom simulations of lOser et alj (|2010L l2012h . These 
simulations were carried out with the following cosmological 
model parameters (in standard notation): h — 0.72, At 
(M'j I t. Adm = 0.216, as = 0.77, and initial slope of power 
spectrum n s = 0.95. First, dark matter-only initial condi- 
tions were evolved from z — 43 to z — 0, and selected indi- 
vidual halos were identified together with their virial radii 
r V ir = J"2oo- F° r this simulation, the softening radius was 
2.52 kpc h~ l . Then these halos were traced back in time, 
and were replaced with high resolution gas and dark matter 
particles. The new halos were evolved again from z = 43 to 
the present day including prescriptions for star formation, 
supernova feedback, gas cooling and a redshift dependent 
UV background radiation. 

From these high-resolution simulations, we here select 
42 galaxies at z=0 which do not have massive satellites 
at this time. The selection is based on the circular veloc- 
ity curves of these systems within 5 effective radii (5 R e ); 
we require that the estimated fluctuations in the CVC in- 
duced by satellites are smaller than 2%. The final simu- 
lated galaxies have stellar masses ranging from 2.7 x 10 10 Mq 
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to 4.8 x 10 11 M Q within 10% r vil . Their typical effective 
radii (defined as the half-light radius within 1 0% r v j r in 
the 2 - dimensional projection) are R c ~ 2-6 kpc (Os er et al.l 
2010). For comparison, the co-moving softening lengths for 
stars and dark matter particles in the resimulations are 
0.4 kpc hT 1 and 0.89 kpc h~ , so the stellar kinematics 
are resolved for R > 0.5 R e in the smallest systems, and for 
R > 0.2 R c in the large galaxies. 

The star for mation model used for the simulations pre- 
sented here (see lOser et ~al\ l20ld for all details) favours 
efficient star formation at high redshift leading to pref- 
erentially spheroidal systems with old stellar populations. 
The simulations do not produce supernova driven winds 
and a model for feedback from central AGN is not in- 
cluded. Therefore, the fraction of available baryons (in 
every halo) converted into stars of the central galaxies 
in the simulated mass range is typically two to three 
times higher than estimates from models matching observed 
galax y mass functions to simulated halo mass functions 
~ 20ld ; iMoster et all l20ld ; iBehroozi etafl 
Possible physical processes respon- 



Guo et al 



ang 



et al.l 



2012) 



sible for this dis crepancy are strong wind-driving feedback 
from SNII (e.g., iPekel fc Siiklll986l; lOppenheimer fc Pavel 
120081 : iGovernato et alj|20ld : iDalla Vecchia fc Schavell2012l ) 
and/or feedback from super-mass i ve b l ack holes (e.g.. 



Croton et ail l200d: 



20ld ; iTevssier et al.ll20lfl ) 



Pi Matteo et all |200S|; iMcCarthv et al) 



Other simulations with similar specifications (weak 
supernova feedback and no AGN feedback) result in 
galaxies with photometric and kinemat ic properties sim- 
ilar to present day e ll iptical galaxies (INaab et al.l 120071 ; 



||n1 

Johansson et al. 200S ; Naab et alj 120091 ; iFeldmann et alj 
20ld ; iJohansson et al.ll2012f ). The galaxies used here are in 
agreement with early-type scaling relations of mass with 
radius and stellar velocity dispersion. In addition, they 
have close to isothermal total mass distributions, simi- 
lar to some observed ellipticals (e.g. [Gerhard et al. 200 ll; 
Koopmans et al.ll200d ; Ich urazov et alj|20ld ;" iBarnabe et al.l 



201 if ), and their observed size evolution between z ~ 2 and 
z = is in agree ment with recent observational estimates 
([Oser et alj|2012h . 

In contrast to the central regions of the simulated galax- 
ies which are well resolved, in their outer regions (R > 2 R c ) 
the particle noise can be substantial. In addition, although 
we have already removed systems having massive satellites, 
many smaller substructures are still present in the outer 
parts of the remaining galaxies. We have decided to smooth 
out these small substructures rather than taking them out 
one by one. However, we have tested in a few cases that 
the results with both approaches are very similar. To re- 
duce fluctuations in the final velocity fields caused by either 
particle noise or small satellites, we use an N-body code 
(NMAGIC, implemented by de Lorenzi et al. 2007) to tem- 
porally smooth the system while integrating the orbits of 
the particles in the gravitational potential, for one circular 
orbit period at 10 R c (~ 10% r v i T ). Since the total mass in 
gas within 10% r v i r in these systems is small (it is one order 
of magnitude smaller than the mass of stars within 10% r v i r 
for all 42 galaxies), we simply fix the gas particles at their 
initial positions while integrating the stellar and dark matter 
particle orbits. 

The NMAGIC code is a spherical harmonics code with a 



made-to-measure (M2 M) algorithm jSver fc Tremainelll99d ; 
Ide Lorenzi et al]|2007f ). Here we use it as a normal N-body 
code. The Poisson solver adopted in NMAGIC is a spher- 
ical harmonic expansion potential solver. It uses a radial 
grid binned logarithmically within maximum radius r max = 
500 kpc, and there are 400 radial bins. The spherical har- 
monics expansion is carried to ? max = 16. This does not 
resolve the small satellites, which are therefore conveniently 
smoothed out during the integration. 

The potential and gravitational acceleration are com- 
puted at a sequence of time steps separated by 



dts : 



1 2vr x 10 R c 



1000 fcirc.10 



R ( 



(i) 



i.e., 1000 times in one circular orbit period for r = 10 R c . 
Puring each interval dts the orbits of the particles are inte- 
grated with an adaptive leap-frog scheme. Kinematic line-of- 
sight (LOS) observables are projected on a two-dimensional 
polar grid with resolution n r x — 20 x 30 on a (10 R c ) 2 
region, and then these observables, A, for the jth grid cell , 
are time averaged by integrating ([Sver fc Tremainel Il99d : 
Ide Lorenzi et al)l2007l ) 



A j (t)=a 



A,(t- 



r)e aT dr, 



(2) 



where the temporal smoothing parameter a is taken to be 



(3) 



1000 dts' 



so that the smoothing time is a" 1 = 200 dts. This procedure 
effectively smoothes over the particle noise for the observ- 
ables in the sparse outer regions of the simulated galaxies. 
In a later section (@, we will illustrate the effect of the 
particle noise and of the satellites on the projected kinemat- 
ics for one snapshot and compare with the time averaged 
kinematics obtained by the procedure just described. 



3 MASS DISTRIBUTIONS AND CIRCULAR 
VELOCITY CURVES 

In this section, we investigate the stellar and dark mat- 
ter density profiles of the simulated galaxies at z=0. We use 
the amplitude and slope of the circular velocity curve to 
characterize the mass distribution, and determine the frac- 
tion of dark matter at intermediate radii. 

Figure [1] shows the three-dimensional volume density 
profiles for stars and dark matter, for radii greater than 
the respective softening radii for the star and dark matter 
particles. The densities shown are temporally smoothed as 
described above. We find for all model galaxies that the 
stellar components have steeply decreasing density profiles 
from 1 R to 6R C , with small cores in the centres, while the 
dark matter halos have flatter density profiles. Since the high 
resolution re- simulations of individual halos include a variety 
of physical processes for the baryonic component, the dark 
matter dens ity profiles are different fr om the simple NFW- 
like profiles (|Navarro et alJll996Ll201ch found in simulations 
that only include dark matter. 
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Figure 1. Temporally smoothed mass density profiles for stellar 
components (magenta lines) and dark matter halos (black) of the 
simulated galaxies. The profiles are shown as a function of nor- 
malized r/R e , extending from the softening length (0.4h kpc 
for stars and 0.89/i — 1 kpc for dark matter particles) to 6R C . 




3.1 Stellar density profiles 

The surface density profiles of observed ETGs can be we ll 
fitted with a Sersic profile l|Sersid Il963l : ICapacciolil ll989T ) . 
Luminous ETGs have 'cored' profiles for which the slope of 
the central profile i s below that of the outer Sercic profile 
( Truiillo et alll2004h We therefore use a cored Sercic profile 
( Graham et al.l "l2003) to represent the stellar density of the 
model galaxies, using 



I(R) 



R' 



7 exp 



l/n 



- 1 



R? 



R 0> 



b n = 1.9992n - 0.3271. 



(4) 



(5) 



Here n is the Sersic index, 7o is a normalization parameter 
related to the central surface density, and R' e is a radius close 
to the effective radius if Ro/R' e <C 1 (because the relation 
used for b n (n) is for the original Sercic profile). To determine 
the parameters of the cored-Sersic profile (Jo, n, Ro and 
R' e ), we fit the density profile out to a truncation radius of 
10% r v i r ~ 10 R c . The particles are binned on a logarithmic 
radial grid (see [J2]), and the surface densities are computed 
and temporally smoothed before fitting the model. 

The cored-Sersic model generally fits the particle dis- 
tributions well, with a typical residual of < 10% of the local 
surface density. Values of the fitted parameters are given in 
Table [2] (columns 2-5) for all model galaxies (designation in 
column 1). We also list the proper effective (2D half-light) 
radius R c for all stella r particles within 10% r v i r (column 
6), which we take from lOser et al.l (|2010l ). We find that R c 
is generally smaller than R' e , with differences < 20% ex- 
cept for the two largest systems. Hereafter, we use R e as 
length scaling unit. Fig. [2] shows the surface density of stel- 
lar particles (black lines) and their cored-Sersic fits (magenta 



Figure 2. Temporally smoothed surface density profiles for the 
stellar components of the simulated galaxies (black lines), to- 
gether with cored Sersic fits (magenta lines), in bins of circu- 
lar velocity at 5 Re. Upper panel: f c i rc (5 Re) < 220 kms -1 , mid- 
dle panel: 220kms~ 1 < i> c ; rc (5Re) < 300 km s -1 , lower panel: 
fcirc(5R c ) > 300 kms" 1 ). 
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Figure 3. Relation between Sersic index and stellar mass for the 
simulated galaxies (left panel) , and relation between Sersic index 
and effective radius (right panel). 
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Figure 4. Left panels: Contribution to circular velocity curves (CVCs) from stars (magenta) and dark matter (black) for all model 
galaxies, separated in bins of increasing total circular velocity at 5R e , as given on the corresponding right panels. The cyan lines show 
the parameterised fits to the dark matter CVCs, using eq. [6] Massive halos with high u c i rc have increasing dark matter CVCs (bottom 
panel), while low-mass halos have flat CVCs (top). Right panels: Total CVCs for all model galaxies: massive galaxies have flat total 
CVCs at large radii (bottom panel), while less massive galaxies (top) have mildly decreasing CVCs. 



lines) . The model galaxies in this figure are binned into three 
groups by mass, i.e., circular velocity at large radii (5R e ): 
fcirc,5Ro = [160,220), [220,300), [300,400) kms" 1 ; see 
below. 

Interestingly, we find from Table [2] that model galax- 
ies with larger stellar masses and effective radii have larger 
Sersic indices, around 6 — 13, while smaller galaxies have 
smaller Sersic indices, around 3 — 6. The left panel of Fig. 
[3] shows the relation between Sersic index and stellar mass 
(within 10% r v i r ) for all model galaxies. We indeed find a 
(weak) correlation between n and stellar mass as described. 
The right panel of Fig. [3] shows the relation between Sersic 
index and size (i.e., R e ). Again we find a weak correlation, 
mostly for the smaller model galaxies; there is a substan- 
tial spread for the larger galaxies. This also implies that the 
sizes and stellar masses of the model galaxies are correlated. 
Both diagrams are also consistent with a bi-modal distribu- 
tion; i.e., a different range of n for systems with effective 
radii R c < 3 kpc and R e > 3 kpc. In any case, systems with 



large R G or stellar mass have on average shallower outer 
density slopes, i.e., more extended outer stellar halos. These 
relations between the Sersic indices, stellar masses and sizes 
of the simulated galaxies are likely to be a consequence of the 
accretion-dominated (major an d minor mergers) l ate assem- 
bly h istory of massive galaxies l|Oser et al.ll2012l ; iHilz et alj 
2012). They ap pear qualitatively consistent with observa- 
tions of ETGs llCaon et al. | 1 19931 : iPrugnieT fc Simienl [l997l ; 



iKormendv eTa l. 2009). 



3.2 Dark matter mass distributions and circular 
velocity curves 

The density profiles for the simulated dark matter halos, 
shown in Figure[TJ are very similar to each other, with nearly 
the same slope between 1 R e and 6R C . To describe the mass 
distributions of the simulated galaxies more quantitatively, 
we study their circular velocity curves (CVCs) in Fig. [4] 
Here the circular velocity v c i IC (r) = [GM^/r] 1 ^ 2 serves as 
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a measure of the spherical part of the mass distribution, in- 
dependent of the actual rotational velocities of the stellar 
or dark matter particles. The three rows of Fig. [3] show the 
model galaxy CVCs binned into three groups according to 
the value of v c i IC at 5 R G (for the binning the same ranges in 
u c i r c are used as in Fig. ^ . The left panels show the CVCs 
for the stellar components (magenta) and dark matter ha- 
los (black) separately. We see that the dark matter CVCs 
are nearly flat at large radii for the lowest v c i rc bin, but in- 
crease outwards for higher mass systems (with high v c i IC ). 
The figure also shows that the dark matter halos dominate 
the CVCs outside 4R C . 

The dark matter CVCs can be well approximated by 
power law mass distributions (potentials) within 10 R e . We 
use the following parametrisation: 



f'n 



(5Rc)° % AT+7 7 ' 



(6) 



Here vo is a normalization parameter for the circular veloc- 
ity, r c is the core radius of the dark halo, and a is the slope 
of the CVC at large radii. Columns 7-9 in Table [T] give the 
best-fitting values of vo, a and r c for the model dark matter 
halos. The values of r c range between 1 — 2kpc for almost 
all models, less than about twice the softening radius of the 
dark matter particles, r^ M = 0.89 ft _1 kpc, and the r c do 
not correlate with the masses of the model galaxies. There- 
fore the presence of a core could be an effect of the softening 
in the simulations. In order to check this further, we also re- 
did the fits while fixing r c = r? M . The best-fitting CVCs 
obtained in this case are not significantly different from the 
previous CVCs (the typical residual is less than 1%). Con- 
sidering the fact that r c /(5R c ) <C 1 (Table [2j, we see that 
the core radii do not significantly affect the rotation curves 
at 5R e . Thus no and a are approximately equal to the value 
and logarithmic slope of the dark matter CVC at 5R C , 



Vo 



!(5Re), 



nDM _ 
Os = 



— 71 (5 R e ). 

a m r 



(7) 



The cyan lines in the left panels of Fig. [4] show the parame- 
terised CVCs of the dark halo components. They agree well 
with the binned data for the halo circular velocities, espe- 
cially at large radii (black lines). 

The left panel of Fig. [5] quantifies the correlation be- 
tween the outer slope S^ M of the dark matter CVC with 
the amplitude of the dark matter CVC at 5 R c for the sim- 
ulated galaxies. We find that for the low-u^^(5R c ) systems 
the slopes of dark matter CVCs are around zero and for the 
high- (5 R c ) systems these slopes are slightly above zero 
(< 0.3), confirming the result from Fig. [4] 

The right panels of Fig. |4] show the total CVCs of the 
model galaxies, including the contribution from the stars, 
which we parametrize again by their value and logarithmic 
slope at 5 R c , 



Vcirc (5 Re 



5, 



d In Weir 

din r 



(5R e ). 



(8) 



The total CVCs are slightly falling at large radii (5 R c ) 
for the systems in the upper right panel, whose u c irc(5Rc) 
is smaller than 220 km s -1 , while they are nearly flat 
outside 2 R e for the most massive model galaxies whose 
Wcirc (5R e ) is larger than 300kms _1 (lower right panel; see 
also lLvskova et al]|2012t ). The right panel of Fig.[5]quantifies 
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Figure 5. Left panel: The logarithmic slope of the dark matter 
CVCs for all model galaxies, versus their outer circular velocity, 
ug^(5Re). Right panel: The logarithmic slope of the total CVCs 
including stars and gas, versus total n c ; rc (5 R e ). Massive galaxies 
(with large values of v c - lrc (5 Re)) have zero slopes, i.e., flat CVCs, 
while smaller galaxies have CVC slopes between flat (0.0) and 
slightly falliing (-0.3). 



this correlation by showing the total slopes 5*5 versus the cir- 
cular velocities at 5R e . Model galaxies with larger values of 
i>circ(5R e ) have S5 ~ 0.0, i.e., flat CVCs at large radii, while 
the remaining galaxies have slopes between [0.0, —0.3], i.e., 
CVCs between flat (0.0) and mildly falling (-0.3). 

Thus it is clear that the slope of the CVC at large radii 
is correlated with the amplitude of the circular velocity, and 
hence with model galaxy mass. The correlation is clearest 
for the dark matter halos alone (see left panel of Fig. [5] and 
Eq. [6}, while it is weakened when the baryonic component 
is taken into account. All values for circular velocities and 
CVC slopes at 5 R e can be found in Table [1] 

The correlation between Sersic index and stellar mass 
shown for the simulated galaxies in Fig. [3] suggests that 
also their total mass and hence circular velocity at large ra- 
dius might correlate with the Sersic index. The upper right 
panel of Fig. [S] confirms such a correlation between the to- 
tal Vcirc (5 R e ) and Sersic index n, but with broad spread. 
Model galaxies with smaller u c irc(5R e ) are on average more 
compact (have smaller n), while model galaxies with larger 
Vcirc (5 R e ) are more extended (have higher n). 

Furthermore, since the circular velocity at large radius 
correlates with its own slope, there should be a correlation 
between the slope of the CVC and the Sersic index. We see 
this correlation in the lower right panel of Fig. [6] Sersic n 
increases with increasing slope of the CVC at large radii. 

At large radii most simulated galaxies are dominated by 
dark matter; therefore, we expect that the circular velocities 
contributed by the dark matter halos also show these cor- 
relations. We show the relation between the Sersic indices, 
Uch^ a t 5R C , and the slope of Uc?rc(5 Rc) in the left panels of 
Fig.[6]and find that these correlations also exist. In fact, they 
are slightly stronger than those obtained with the parame- 
ters of the total mass distribution which includes the stellar 
and remaining gaseous mass in the simulated galaxies. 

Finally, we show the fraction of dark matter within 3- 
dimensional spheres of 1 R c (left panel of Fig. [7} and 5 R c 
(right panel) for the model galaxy sample. At 1 R G , the dark 
matter fractions are between 15% — 30%, i.e., the luminous 
matter dominates by a factor 3-7. At 5R C , the dark mat- 
ter fractions are between 40% — 65%, i.e., the amount of 



© 0000 RAS, MNRAS 000, 000-000 



The mass and angular momentum distribution of simulated massive galaxies to large radii 7 



15 



10 



A A A 

a m .a 



*a. a a A . 



V A^ A 




100 150 200 250 300 550 150 200 250 300 350 400 
v c:rc.5Re DM (km/s) v cirCi5Re (km/s) 




-0.2 -0.1 0.0 0.1 0.2 0.3-0.3 -0.2 -0.1 0.0 0.1 
Slope of v drCi5Re DU Slope of v clrc5Re 

Figure 6. Relation between Sersic index and circular velocity 
curve parameters for simulated galaxies. Sersic index is plotted 
against circular velocity from dark matter only (upper left panel), 
slope of the dark matter CVC (lower left panel), total circular 
velocity (upper right panel), and slope of total CVC curve (lower 
right panel) , all evaluated at 5 R c • 



dark matter is now on average larger than the mass in stars. 
We can infer from Fig. [7] that the DM fraction increases 
with stellar mass (or circular velocity, or Sersic n), and that 
this increase is mostly due to a dependen ce on effective ra - 
dius rather than circular velocity (see also lHilz et al.|[2012t ). 
These results are qualitatively consistent with dark matter 
fractions inferred from po wer-law dynamical m odelling of 
PNe an d GC observations l|Deason et al.l [20121 ) and strong 
lensing (jAueer et al.ll2010l ). 

3.3 Shape correlation and alignment 

Observed ETGs have ellipticities up to e ~ 0.8 
(jBernardi et al.ll2003l : Ikra movie et al.ll201ll ). The distribu- 
tion of ellipticities depends on sample selection; it is approx- 
imately flat up to ellipiticity e ~ .7 for the recent Atlas3D 
sample (see iKrainovic et al 1 l201ll . Fig.7). SO- like fast rota- 
tor ETGs are consistent with being a family of oblate sys- 
tems viewed at random inclination angles, while slow rota- 
tor ETGs are moderately flattened (e ~ 0.3) and often show 
kinemati cal misalignments, ind i cating that they are tr iaxial 
systems jCappellari et alJ^OoilEmsellem et al.ll201ll ). 

The shapes of the outer dark halos (> 10% r v i r ) can 
be estimated from the shear patterns in weak gravitational 
lensing dat a (Hoekstra et al. 2004; Mandelb aum et al]|2006l ; 
van Uitert et al.ll2012l ). Some of these studies have found 
that the dark matter halos of red galaxies on scales beyond 
~ 0.1r v i r are preferentially aligned with the lens galaxies, 
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Figure 7. Dark matter fraction for all model galaxies at 1 R G (up- 
per panels) and 5R C (lower panels), versus outer circular velocity 
' u circ(5Re) (left panels) and versus effective radius (right panels). 
Within lR e , the average mass in dark matter is about 25% of 
the total, while within 5R e , the average dark matter fraction is 
slightly above 50%. The dark matter fraction has little (at 5 R c ) to 
no dependence (at 1 R e ) on circular velocity, but increases clearly 
with R c . 



but the signal is not as clear as might have been expected. 
Dark matter simulations predict that t he halos are triaxial 
(|jing fc Sutoll2002l ; lAllgood et al.ll2006l ). In hydrodynamical 
simulations of disk galaxies, the inner halos become prefer- 
ent iall y aligned with the disk and misaligned with the oute r 
halos (Bailin et al.ll2005l ; iBett et al.ll201ol ; lHahn et af1l2010l ), 
while in binary major mergers the short axis of the stellar 
remnant is found to be orie nted perpendicular to the long 
axis of the surrounding halo (|Novak et al.ll2006l ) . It is there- 
fore of interest to investigate the same question for the cur- 
rent model galaxies, which have a very different formation 
history. 

To study the alignment of the simulated galaxies and 
their host dark matter halos, we determine the principal axis 
directions of both components at small radii (2 R c ) and large 
radii (5 R c ) at z~0. We find these quantities by diagonalising 
the moment of inertia tensor iteratively from all particles in- 
side ellipsoids x 2 + (y/p) 2 + (z/q) 2 = f 2 R e 2 , where / = 1 or 
5, p — b/a, q = c/a and, e.g., the latter is determined from 



q = c/a = ulzz/lxx where I xx , hz are the diagonalised 
moments of inertia for all stellar particles within the ellip- 
soid. Because within 1 R c there are only a few hundred dark 
matter particles for many of our models, we have chosen 
the inner radius at 2 R c in order to ensure reliable results. 
The particle mass for the stellar component is ~ 1/5 of the 
dark matter particle mass, so there is a sufficient number of 
stellar particles even within quite small radii to define the 
shapes of the 'luminous' galaxies. 

Figure [8] shows the histograms of misalignment angles 
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Figure 8. Alignment of the principal axes of the stellar and dark 
matter components. Histogram of misalignment angles for all 42 
simulated galaxies, between: principal axes of the stellar compo- 
nent at 2 R and 5 R e (top left) ; of the dark matter component at 
2R C and 5R C (top right); of stars vs. dark matter at 2R C (bottom 
left); and stars vs. dark matter at 5R C (bottom right). Colours 
in each panel denote misalignment angle of the respective short 
axes (cyan), intermediate axes (magenta), and long axes (black). 



between the short, intermediate, and long axes of the stellar 
and halo components at 2 R c and 5 R e for each component 
separately, and between both components. The short axes 
(blue) of the stars at different radii, the dark matter at dif- 
ferent radii, and of the stars and dark matter distributions 
are always very well aligned (within < 5°). The long axes 
(black) and intermediate axes (purple) are approximately 
aligned within (< 15°) in about half of the cases; in the 
other half of the cases they are misaligned by roughly 90 
degrees, meaning that the intermediate and long axes have 
switched between different radii or components. This can 
happen most easily when the system is nearly axisymmet- 
ric. 

Figure [9] shows the short-to-long axis ratios of the dark 
matter halos versus those of the stars, at both 2 R c and 
5 R c . These axis ratios where obtained by viewing the system 
along the intermediate axis of the stellar component, and 
computing the projected axis ratio c/a from the diagonalized 
2D projected moment of inertia tensor. At both 2 R c and 
5 Rc, the axis-ratios of the galaxies are in the range [0.4,0.8] 
while the halos are slightly rounder, c/a £ [0.5,0.95]. From 
the figure we see a moderately strong correlation between 
the shapes of the stellar and halo components, with scatter 
in the halo axis ratio at given stellar c/o of ~ ±0.1 at 5R C 
and ~ ±0.2 at 2R C . 

In summary, the short axes of the simulated galaxies 
and their host dark matter halos are well aligned within 
< 5° througout the radial range probed (2R C -5R C ), and 
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Figure 9. Short-to-long axis-ratios of stellar components versus 
halo components at 1 R e (left) and 5R C (right). The c/a of stars 
and halos strongly correlate at large radii but not at small radii. 



their shapes are correlated. Long and intermediate axes are 
either aligned or misaligned by 90°, i.e., switch their order, 
within slightly larger (< 15°) scatter. 



4 LINE-OF-SIGHT KINEMATICS OF 

SIMULATED GALAXIES TO LARGE RADII 

As is well-known, the outer kinematics of ETGs are difficult 
to measure because the stellar densities and surface bright- 
ness profiles decrease rapidly at large radii. The same prob- 
lem also exists in the simulated galaxies where the particles 
follow similar density profiles. In addition to Poisson noise, 
there is also a further source of fluctuations that arises from 
various small satellites around the central galaxies. 

In order to investigate the effects of the fluctuations 
from both low particle numbers and satellites, we consider 
three cosmological galaxies from our sample of lOser et alj 
(|2010h more closely. These are denoted M0125, M1017, and 
M0300, and are selected based on their rotational proper- 
ties (see also § [5} ■ Velocity maps for these galaxies at z = 
are shown in Fig. 1101 The projection direction is along the 
intermediate axis of the stellar distribution within one effec- 
tive radiufl The unsmoothed snapshot maps of mean LOS 
velocity v (upper panels in the left column of Fig. 1 10 p show 
significant fluctuations. A map of the error 5v in the mean 
LOS velocity is shown in the lower left panels of Fig. 1101 
This is defined as Svj = aj/y^Nj, where <jj and Nj are the 
velocity dispersion and number of particles in the j t h cell 
grid for this line-of-sight. Typical values for these fluctua- 
tions are 5v 3 = 20-80 kms" 1 . Model M1017 is an extreme 
case where the error can reach almost 50% of the mean LOS 
velocity itself. 

In addition to the fluctuations caused by particle noise, 
these velocity maps show a number of well-defined substruc- 
tures. These satellites will be moving along their own orbits 
until they finally merge with the central galaxy, and can lo- 
cally have quite different LOS velocities from the host galaxy 
particles. Thereby the satellites change the 1-dimensional 
LOS velocity profiles and thus v, and they may also affect 
the angular momentum profiles. We have decided not to 



1 The maps are interpolated from a polar grid as described in 
Sect. [2] 
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Figure 10. Edgc-on LOS kinematics of selected model galaxies. Upper panels, left column: mean LOS velocity maps for the stellar 
compone nts of three cosm ological galaxies (labeled M0125, M1017, M0300), computed from the particle distributions for the snapshot 
at 2 = llOser et al.ll201(t) . The maps are interpolated from a grid in radius and angle of n r X = 20 X 30 cells on a (10 R e ) 2 region. 
Signatures from both particle noise and from a number of small substructures are visible. Lower panels, left column: particle noise error on 
the mean velocity for the same galaxies. The grid points (black dots) are overplotted. Upper panels, right column: Temporally smoothed 
mean LOS velocity maps (v) for the stellar components of the same model galaxies. Middle panels, right column: Maps of temporally 
smoothed velocity dispersion rr. Bottom panels, right column: Maps of ratio v/a. The surface densities of the stellar components are 
overplotted on these maps; contour levels are 10 70 , 10 7 ' 5 , 10 80 , 10 8 ' 5 , 10 9 , 1O 9 ' 5 M0 kpc 2 . The three model galaxies chosen here are: 
a slow rotator (M0125), a slow rotator with a peak of X(R) around 2R C (M1017), and a fast rotator with increasing A(i?) up to 2R e 
(M0300); see Sect. [5] 



take all these small satellites out, but to consider the host 
galaxy after these substructures have phase-mixed away in 
its large-scale gravitational potential. The NMAGIC code 
used for this purpose (see Sect. [2j also computes a time 
average of the observables v and Oj which allows us to ob- 
tain more reliable LOS kinematics at large radii than would 
otherwise be possible. 

Recomputing in this way the projected LOS kinematics 
for the three galaxies in Fig.[Trj]results in the kinematic maps 
presented in the right column of Fig. 1101 The three maps in 
the upper right panel are mean velocity maps v. Fluctuations 
are greatly reduced compared to the unsmoothed maps in 
the left column. Also the fluctuations arising from the small 
satellites are mostly smoothed out. The middle panel in the 
right column of Fig. [TU] shows maps of LOS velocity disper- 
sion a, and the lower panel shows maps of the ratio v/a with 
overplotted stellar surface density contours. 

For model M0125 (always in the left panels), the ve- 
locity field shows a mild rotation, and its velocity disper- 
sion map shows large a within 2.5 R e , implying that this 



galaxy is pressure-supported. For model M0300 (right pan- 
els), rapid disk- like rotation is seen from the mean velocity 
map, while a is only large in the centre (R < 0.5 R c ); this ob- 
ject is significantly supported by rotation. Finally, in model 
M1017 (middle panels), one can see a hot disk- like structure 
at intermediate radii (~ 2.5 R c ) in both the velocity and 
a maps, while the rest of the galaxy has lower dispersion. 
Model M1017 is mostly pressure-supported with a rotation- 
ally supported component. While this model is somewhat 
unusual, the other two are quite typical for the large radius 
kinematics of fast and slow rotator galaxies in the simulated 
ETG sample. Observed outer velocity fields similar to the 
three cases shown here are those of NGC 5846 (a slow ro- 
tator galaxy), NGC 4564 (a fast rotator with a disk-like ve- 
locity field at intermediate radii; for both see lCoccato et all 
2009), and NGC 1316 (a rapid ly rotating merger remnant; 
see lMcNeil-Movlan et ai1l2012l ). 

Figure [TT] shows line-of-sight velocity fields for 6 ad- 
ditional simulated galaxies out to 5 R c . These galaxies have 
been chosen as representatives for different dissipational and 
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Figure 11. Edge-on mean LOS velocity maps for the stellar com- 
ponents of further 6 selected model galaxies. Upper panels show 
three simulated galaxies with large scale rotation. The galaxy in 
the top right panel has a pronounced disk component from a 
gas rich major merger. The lower panels show three simulated 
galaxies with a dissipationless merging history. The galaxy in the 
bottom right panel shows weak major and minor axis rotation; it 
has a formation history with many minor mergers. 



dissipationless formation histories. For example, the galaxy 
in the top right panel has a pronounced disk component 
around R ~ 2R C from a gas rich major merger, and the 
galaxy in the bottom right panel shows weak major and 
minor axis rotation; it has a formation history with many 
minor mergers. For further details about the formation his- 
tories of these simulated galaxies see Naab et al. (2012, in 
preparation) . 



4.1 Root Mean Square Velocity Profiles 

Observations of outer ETC kinematics using planetary neb- 
ula velocities have shown that most ETGs are characterized 
by slowly declining profiles of circularly averaged RMS LOS 
velocity v IYns (R), but with a significant mino rity of galaxies 
for w hich the RMS velocity declines rapidly l|Coccato et alj 
2009). It is therefore of interest to investigate the equivalent 
radial profiles for the (re)simulated galaxies considered here. 

Figure [T5] shows the circul arly averag e d pro files of 
t) rms (fl) for the 42 galaxies from lOser et al.l (|2010T ) in the 
three usual bins of u c irc(5R c ). The viewing direction for all 
simulated galaxies is edge-on along the intermediate axis of 
the stellar distribution. The w rms (rj) include projected rota- 
tion and velocity dispersion of all star particles, and are nor- 
malised by the respective values of u rms (lR c ). We find that 
the v Ilna (R) profiles decline moderately with radius for the 
low-mass systems, whereas they decrease only mildly for the 
high-mass systems. This is consistent with the CVCs shown 
in Fig. [4] The v Yln a(R) profiles of fast and slow rotators (see 
Sect. [5]) are not significantly different. 
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Figure 12. The root mean square velocities » Ims of the stellar 
components of 42 galaxies versus radius in bins of u c i rc 5 ri e • The 
t) rms (fl) profiles are normalized by their respective v TulB at 1R C . 
Solid lines are for fast rotators while dotted lines are for the slow 
rotators; see Section [5] Coloured curves depict the three galax- 
ies with velocity fields shown in Fig. 1101 M0125 (green dotted 
line, bottom panel), M1017 (blue dotted line, middle panel), and 
M0300 (red solid line, middle panel). 



The Vrms( R) profiles in Fig. | 12l are consistent with the 
major group in lCoccato et al.l l|2009t) . but there is no equiva- 
lent in the simulations for the rapidly falling profiles shown 
in that paper (e.g., for NGC 3379). For this comparison, 
we also studied the u rms (R) profiles of the simulated galax- 
ies in face-on projection, with LOS parallel to the shortest 
axis. These face-on profiles are similar to those in edge-on 
projection. At least for the current sample of cosmological 
galaxies, the inclination angle does not seem to be an im- 
portant parameter for n rms (i?). However, following the for- 
mation of thin disks is still problematic for such simulations, 
and this may also impact the geometric properties of merger 
remnants from such disks. Therefore this question should be 
revisited with future simulations once the disk problem has 
been solved. 



4.2 Anisotropy profiles 

In these cosmological simulations, the stars in the model 
galaxies have two origins: the inner in situ component forms 
at early times in a period of rapid star formation, while 
the second component stems from accreted stellar particles 
which is predomi nantly but not only found at large radii 
(|Oser et al.l I20l6l ). The accreted component is character- 
ized by radially anisotr opic velocity dispersions l|Abadi et alj 
120061 : lmiz et al.ll2012T) because the merging satellites come 
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Figure 13. Velocity anisotropy profiles for the simulated galaxies 
in bins of in situ fraction of stars within 10%r v i r . The coloured 
curves are as described in Fig. 1121 Galaxies with low in situ frac- 
tions have more radially anisotropic orbit distributions, while tan- 
gential anisotropy is seen only for systems with high in situ frac- 
tion. 



in on predominantly radial orbits, and so many of the stars 
stripped from the satellites enter the host galaxy on nearly 
radial orbits also. 

Therefore, we would expect the kinematics of model 
galaxies to be more or less radially anisotropic depending 
on the relative fractions of in situ and accreted stars. Figure 
1 13 1 shows anisotropy profiles f3(r) si — (ag + a^)/2a^ for all 
model galaxies, divided into bins of in situ fraction. As ex- 
pected, tangentially anisotropic model galaxies can only be 
found in the group with large fraction of in situ stars (lower 
panel and also M1017 in the middle panel). Almost all sim- 
ulated galaxies with low in situ star fraction (upper panel) 
are radially anisotropic with j3 ~ 0-0.3; but many of the sys- 
tems with higher in situ fraction have similar anisotropies. 
As a result, there is no correlation between both quantities 
for the whole sample. We also find that slow rotators are 
radially anisotropic (except for M1017), while fast rotator 
model galaxies have both radially and tangentially biassed 
anisotropy profiles. 



5 ANGULAR MOMENTUM PROFILES 

To characterize the s pecific angular momentum of ETGs, 
lEmsellem et al.1 (^OOTT l introduced the X{R) profile X(R) =< 
Rv > I < R\/v 2 + a 2 >, which is a luminosity- weighted, 
cumulative measure of projected angular momentum per 
unit mass within radius R. They found that ETGs can be 



divided into fast and slow rotators, according to whether 
X(R = R e ) > 0.1 resp. < 0.1. Analyzing simulated merger 
remnants, Ijesseit et al.l (|2009l ) found that \(R) is a good 
proxy for the t r ue an gular momentum of these remnants. 
lEmsellem et al.1 (|201lh revised the criterion for separating 
slow and fast rotators to include ellipticity values at R e or 
R E /2, and they also considered the effect of inclination and 
an assumed anisotropy on the resulting X(R) profile. 

Much of the angular momentum of ETGs could reside at 
large radii, where kinemati c measurement s are m ore difficult 
than in the bright centers. ICoccato et al.l (|2009l ) studied the 
kinematics of the outer halos of a sample of ETGs with PNe. 
They found (their Figure 14) that the X(R) profiles to large 
radii for the most part confirm the separation into fast and 
slow rotators based on the inner kinematics, but with some 
slow rotators having up to X(R) ~ 0.3 at large R and some 
fast rotators whose X(R) profiles decrease outwards. 

It is therefore of some interest to compare the specific 
angular momenta of the simulated galaxies with those of 
observed ETGs, out to several effective radii. To obtain an 
overview of the angular momentum properties of the sim- 
ulated galaxies, we study the X(R) profiles for all the 42 
cosmological galaxies, and to reach large radii we use the 
temporally smoothed observables (Sect. [2} in the computa- 
tion of X(R), 
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where the summation is over the kinematic grid. 

Figure [14] shows the resulting X(R) profiles, again with 
the simulated galaxies binned according to their outer cir- 
cular velocity i> c irc,5R . The left panels of Fig. 1141 show the 
X(R) profiles for edge-on view, while the right panels show 
the face-on X(R) profiles. Most of the simulated profiles have 
a very regular form, reaching a nearly constant value at 
~ 1 — 2R e , but some increase or decrease more noticeably 
at large radii. Both these variations and the general shapes 
of the X(R) profiles from the cosmological galaxies appear 
consistent with the observed X(R) p rofiles of ETGs to large 
radii, shown in lCoccato et al.1 (|2009h . Their slopes at small 
radii appear somewhat shallower than the typical slopes seen 
in th e SAURON and ATLAS 30 data l|Emsellem et alj|2007l . 

1ml). 

To obtain a more detailed picture of the angular mo- 
mentum in the outer stellar halos, it is useful to study the 
more local angular momentum parameter \v\/a(R) as well. 
The local \v\/a profiles for all simulated galaxies in edge-on 
view are shown in the middle panels of Fig. 1141 For most 
simulated galaxies, the \v\/a(R) profiles are nearly flat or 
increasing mildly with radius after reaching a plateau at 
1 — 2R e . However, there are also a few exceptions for which 
\v\/a(R) decreases at large R, and some with a strong local 
concentration of specific angular momentum. These features 
do not show up as well in the cumulative X(R) profiles. 

As mentioned above, ETGs can be classified into fast 
and slow rotators according to their X(R) v alues at R e . For 
classi fying the simulated galaxies we follow l|Emsellem et al.l 
120071 ). taking A(R C ) > 0.1 to define fast rotators (solid lines 
in Fig. I14jl and A(R C ) < 0.1 for slow rotators (dotted lines). 
The main reason for this choice is that this classification 
does not depend on ellipticity. We find that among 42 cos- 
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Figure 14. Specific angular momentum parameter A(R) profiles for the sample of 42 cosmological re-simulation galaxies, binned by 
circular velocity f c irc(5R c ) and plotted outside the softening radius OA/h kpc. Edge-on profiles (inclination i = 90°) are shown in the 
left panels and face-on profiles (i = 0°) in the right panels. The middle panels show \v\/a profiles for i = 90° in the same bins; these 
provide more local rather than cumulative information on the rotational properties. The solid lines are for fast rotators whereas the 
dotted lines show slow rotators; this classification depends on viewing angle. The coloured curves point to three typical galaxies: M0125 
(green dotted lines on bottom panels), M1017 (blue dotted lines on middle panels) and M0300 (red solid lines on middle panels). 



mological galaxies, 23 are fast rotators (based on the edge-on 
profiles) and 19 are slow rotators. The fraction of slow ro- 
tators is largest in the group of simulated galaxies with the 
highest circular velocities at 5R e . In the face-on view (right 
panels of Fig.ll4[). none of these systems rotates rapidly, and 
they are all classified as slow rotators. The fraction of slow 
rotators among the simulated galaxies ( 45%) is significantl y 
more than in the ATLAS 3d sample (jEmsellem et , all 1201 fl ); 
how ever, the sample selection criteria are not comparable 
(see lOser et~ai1l2010l ). 

In Fig. [TO] we showed the velocity fields of three typical 
galaxies; a slow rotator with very little rotation at any radii 
(M0125, left panels), a fast rotator for which rotation is 
significant at all radii (M0300, right panels), and a system 
rotating slowly in the inner R e but rotating rapidly for R ~ 
2 — 4R C (M1017, middle panels). Especially the v/a maps 
in Fig. [10] show the different rotational properties of these 
model galaxies well. 

We use the same three simulated galaxies to show the 
beneficial effect of the temporal smoothing on the cumula- 
tive X(R) and local \v(R)\/a profiles in Fig. 1151 where the 
smoothed profiles from Fig. [14] are shown in red and the 
unsmoothed profiles are overplotted in black, for compari- 
son. The upper panel shows that the amplitudes of the un- 
smoothed \(R) profiles are significantly higher than those 
of the time-smoothed counterparts. The lower panels show 
a similar effect in the local |«(-R)| and \v\/a(R) profiles; these 
have large fluctuations for a single snapshot, especially for 



the slow rotators, with values of up to a factor of 2 larger 
than in the smoothed profiles. 

The reason for the higher amplitudes in the unsmoothed 
case lies in the definition of the X(R) parameters in terms 
of absolute values of mean velocity, so that negative and 
positive fluctuations cannot cancel in the angular or radial 
summation. The unsmoothed profiles can therefore be bi- 
ased significantly by particle noise, especially for the slow 
rotators (such as model M0125) and in the center, and by 
global asymmetries (such as seen in model M0300 out to 
large radii). The light-blue curves in the top panel of Fig. 1141 
show the results of an additional test for this effect. They 
are determined from computing X(R) separately for the pos- 
itive and negative velocity sides of the mean velocity maps in 
Fig. 1101 but without the absolute value convention of eq. [9] 
and then adding the absolute X(R) values from both sides. 
This has the effect that positive and negative velocities on 
each side are allowed to cancel, which mimicks the effect of 
removing the noise or the asymmetry by the phase-mixing 
that occurs during temporal smoothing. It is clear from the 
Fig-HHthat the profiles obtained from this asymmetric aver- 
aging and from temporal smoothing agree closely with each 
other for all three simulated galaxies. This shows that to 
obtain a correct indication of the angular momentum of the 
system, the time averaging is necessary. 

Returning to Fig. 1141 it is clear that most of the sim- 
ulated edge-on X(R) profiles are flat to slightly rising for 
R > 2 R c (Figs. [H [15]). Thus, most fast rotators at R e 
continue to rotate rapidly at large radii. Most of the slow 
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Figure 15. The 1-dimensional profiles of mean absolute veloc- 
ity (lower panel), local rotation parameter \v(R)\/a (middle), 
and cumulative specific angular momentum parameter A(R) (up- 
per panel), for the three simulated galaxies M0125 (solid lines), 
M0300 (dashed lines) and M1017 (dotted lines). Black curves 
show profiles computed from unsmoothed data, magenta curves 
show the results based on temporally smoothed data, and light- 
blue curves show profiles based on two-sided averaging (see text). 
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Figure 16. Relation between rotation, ellipticity and mass for 
outer stellar halos and central regions of the 42 simulated galaxies. 
Upper left panel: outer local \v\/a (5R ) versus inner \v\/a (R e ) 
for edge-on view. Upper right panel: A(5R e ) versus A(R e ). Stel- 
lar halo rotation correlates with central rotation, especially for 
the flattest galaxies. Lower panels: X(R) versus ellipticity e, at 
R e (left) and 5R e (right), for edge-on view. The colour code in 
all panels shows stellar mass enclosed within 10% r v i r . The more 
massive galaxies are rounder and have lower specific angular mo- 
menta, while the less massive galaxies are more evenly distributed 
in both parameters. 



rotators have mildly increasing X(R) in the outer regions, 
especially for the group of less massive systems with circu- 
lar velocities at large radii < 220 km s -1 . Fig. 1161 confirms 
these trends. The upper left panel of Fig. [16] shows a corre- 
lation with much scatter between the central and outer local 
specific angular momentum, \v\/a at 1 R c and 5R C . The up- 
per right panel of Fig. [16] shows the close relation between 
the X{R) values at lR e and 5R C , consistent with the flat 
X(R) profiles. The correlation between the two X(R) values 
is stronger than for the local angular momentum parameters 
\v\/a at the same radii, because of the cumulative nature of 
X(R) within radius R. 

The \{R) profiles for a large number of nearby ETGs 
are sh own out to < 1.5 R e in Figure 5 of lEmsellem et all 
(l201ll). More e xtend ed X(R) profiles f or ETGs are sh o wn in 
iProctor et ail l|2009l . out to 3 R e ) and lCoccato et al.l (|2009l , 
out to 9 R ). Comparing with these observations, the X(R) 
profiles of our model galaxies are consistent with observed 
ETGs. Even for the unusual Mo del M1017, within 1.5 R e , 

there are similar A(i?) profiles in lEmsellem et ah (12011 ). 

The observations of the SAURON (ICappellari et al.l 

120071 ) and ATLAS 3d projects (|Emsellem et al.l l201ll ) have 
shown evidence for correlations between specific angular mo- 
mentum parameters (either \v\/a or X(R) profiles) and el- 
lipticity of these ETGs. To compare with these observations, 
we show in Figure [16] the \(R) parameters versus the ellip- 



ticities e of the 42 simulated galaxies, both at 1 R e (bottom 
left) and 5R e (bottom right). Here ellipticity e = 1 — c/a, 
where a, c are the long and short semi-axes of the projected 
ellipse in edge-on projection (see a]3.3[) . We list the values of 
ellipticity at 1 R e and 5R e in columns 10 and 11 of Table 
[5] recalling that the central ellipticities for the lower mass 
systems have some uncertainty due to softening. Some of the 
galaxies contain highly flattened disk-like structures which 
result in large ellipticities at R ~ R c . 

In general, the simulated slow rotators are rounder than 
the fast rotators, similar as in the observations. The distri- 
butions at 1 R e and 5 R c are similar, despite the already 
mentioned possibility that the ellipticity within R c could 
be somewhat influenced by softening. The bottom left panel 
of Fig. [16] is in general agreement with the ATLAS 3d ob- 
servations; however, the ellipticity range is narrower for the 
simulated galaxies. This is because the ATLAS 3d sample 
contains some very flat ETGs which may be disk galaxies, 
while our sample galaxies are mostly spheroidal galaxies. 

The colour bar in Fig. [16] labels the stellar mass of the 
simulated galaxies within 10% r vlT . This shows the corre- 
lations between angular momentum, ellipticity and stellar 
mass for the simulated galaxies: more massive galaxies rotate 
more slowly (have smaller X(R)) and they are rounder, while 
the specific angular momenta for the less massive galaxies 
are spread over a wider range, in broad agreement with ob- 
servations. 
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6 CONCLUSIONS AND DISCUSSION 

In this work, we analyzed the structure and kinematics of 
4 2 simulated galaxi es from the cosmological re-simulations 
of lOser et all (|201(^ which do not contain massive substruc- 
tures at a level of > 2% perturbation in the CVC at redshift 
z — . Here we summarize and discuss our main results. 

1. The stellar components of the model galaxies can be 
well fitted with cored Sersic profiles. The Sersic indices n 
correlate broadly with the stellar masses and sizes of the 
galaxies, with significant scatter, such that more massive 
galaxies have larger n, similar to observations of ETGs (e.g., 
Caon et al.lll9931 ; iPrugniel fc S imicn 1997; Kormendv et al.l 
20091 ). 

2. Due to the interaction between dark matter and 
baryons during the assembly processes, the dark matter 
density profiles in the model galaxies deviate from the 
NFW-like profiles characte ristic for dark matter only models 
(|Navarro et alJI 19961 l2O10h . The DM density slopes and the 
slopes of the DM CVCs are consistent with power laws and 
vary systematically with mass, such that the DM CVC is ap- 
proximately flat (-R ) for less massive systems and slightly 
rising (R 3 ) for high mass galaxies. 

3. The corresponding total CVCs are slightly falling 
(R~ 3 ) for the less massive systems and approximately flat 
(R°) for the more massive model galaxies. This is consistent 
with mass determinations for local galaxies (see the Intro- 
duction). As is then expected, galaxies with higher Sersic 
indices n are found to have larger values of ti c h- c (5R e ) and 
of the slope of the CVC at 5R C . These outer parts of the 
cosmological galaxies are mostly dominated by their dark 
matter halos; so these trends hold also for the dark matter 
CVCs separately. 

4. The short axes of the simulated galaxies and their 
host dark matter halos are well aligned within < 5° 
througout the radial range probed (2R -5R ), and their 
shapes are correlated. Long and intermediate axes are ei- 
ther aligned or misaligned by 90°, i.e., switch their order, 
within slightly larger (< 15°) scatter. 

5. We computed mean velocity v, velocity dispersion a, 
and local v/a fields out to 5 R c for the simulated galaxies 
and illustrated their kinematic diversity for three of these. 
We temporally smoothed these velocity fields in order to 
suppress particle noise and fluctuations from small satel- 
lites and illustrate the necessity and effect of this in some 
detail. The simulated galaxy sample contains both purely 
dispersion-supported systems with no or little rotation, and 
objects that show disk-like rotation at a level of v/a ~ 1. We 
also show a rarer case with a disk- like structure at > 2.5 R c . 
The observed outer velocity fields of galaxies from PNe show 
a similar kinematic diversity. 

6. Radial profiles of root mean square velocity v TUls (R) 
are slowly declining, independent of whether the simulated 
galaxies are fast or slow rotators, and similar for edge-on 
and face-on projections. These profiles resemble the majority 
group of outer v Ilns (R ) profiles determined in nearby ETGs 
from PNe velocities bv lCoccato et all (|2009h ; however, there 
are no analogues in the simulated galaxies for the rapidly 
falling Urms(-R) profiles seen in the observed sample. This 
could mean that these objects form through a channel dif- 
ferent from the simulated galaxies studied here, or it could 



be related to the difficulty of modelling disk galaxies in the 
cosmological simulations; this issue requires further study. 

7. Simulated galaxies with a large fraction of accreted 
stars are generally radially anisotropic. Only systems with 
dominant in situ fraction show tangential anisotropy. These 
trends are due to the fact that in the simulation the ac- 
creted stars are tidally dissolved from merging satellites on 
preferentially radial orbits. We also find that the slow ro- 
tators amongst the simulated galaxies are mostly radially 
anisotropic, while the fast rotators have both radial and tan- 
gential anisotropy profiles. 

8. We determined cumulative X(R) and local v/a(R) 
angular momentum parameter profiles for the stellar com- 
ponents of all simulated galaxies from the time averaged ve- 
locity fields. For most simulated galaxies, the edge-on A(i?)- 
profiles are flat or slighty rising within 2 R c — 6 R c . Most fast 
rotators rotate fast at both small and large radii, but some 
have decreasing rotation at large radii. Lower mass slow ro- 
tators have mildly increasing \(R) with A(6R e ) G (0.1, 0.3], 
whereas high-mass slow rotators have flat A(i?)-profiles. 
Overall, A increases with ellipticity, but with much scat- 
ter. These properties appear broadly consistent with those 
of observed ETGs. 
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Table 1. Stellar mass and circular velocity curves for the simulated galaxy sample: designations of the model galaxies (1st column), 
stellar mass within 10% r v j r (2nd column), circular velocities and their slopes at radius 5R e (3rd and 4th columns), circular velocities 
from dark halo components only and their corresponding slopes at 5R e (5th and 6th columns), parameters of the dark matter halo fit of 
eq. [6](7t/i to 9th columns). 



Model Mass* ^ c irc,5R c Slope of t)™ 5B Slope of vq a r c 

(W 10 M Q ) (kms- 1 ) » circ , 5Ho (fan's" 1 ) «™ 5Re (kms" 1 ) (kpc) 



Mnno4 
iviuuy4 


A 7 
4 1 . 


on 
.yu 


■388 
• >oo . 


8i 

Ol 


u 


09 
.UZ 


904 

zy4. 


'39 
. OZ 


u 


9 A 
. Z4 


9QA 

zyo. 


no 
. uu 


u. 


90 
. ZU 


1 

1 


89 
.OZ 


1V1U1ZO 


A '3 
4o. 


■3 a 

-OU 


OiO. 


on 
. OO 


u 


O Q 
.Uo 


90S 

Zyo. 


'3 O 
. OU 


u 


9 '3 

.Zo 


907 

zy ( . 


AO 

. uu 


U. 


1 O 

. iy 


1 

1 


Q7 

.y f 


1V1U1 I 


'3 A 


70 
. ( y 


*379 
O / Z. 


A A 
.11 


n 
-u 


0'3 
.UO 


9(H) 

zyy. 


A r ! 
■IO 


n 


1 A 
. 14 


9QA 

zyo. 


1 8. 

. IO 


n 
u. 


1 A 
IO 


o 
Z 


09 
.UZ 


Mm on 
iviuiyu 


■3 "1 

ol . 


A fi 


oZU. 


. lo 


n 

-U 


.Uo 


9 A "1 
Z41 . 


"1 A 
. ID 


n 
u 


OO 

.uy 


9 1 - 3 
Z40. 


A 9 
.OZ 


U- 


1 - 3 

. lo 


1 
1 


9Q 
.Zo 


1V1UZU4 


ZO. 


8 A 
. o I ) 


ouo. 


OO 


u 


'3 
.UO 


99Q 

zzy . 


'3 1 

. Ol 


u 


97 
. Z f 


9*30 
ZOU. 


70 


n 

u. 


99 
. ZZ 


1 

1 


i n 

. 1U 


mooog 
iviuzuy 


1 Q 

iy . 


.yo 


ouy. 


1 A 

. lo 


-U 


O O 

.uy 


9Q9 
ZOZ. 


A 7 
.4 / 


u 


1 '3 

. lo 


9*3A 
ZOO. 


("17 
A) ( 


U. 


OA 
UO 


z 


i 1 

.14 


1V1UZ10 


97 
Z i . 


a-1 


OO ( . 


. OU 


n 
u 


1 
.U 1 


9^0 

zoy 


'39 
. OZ 


n 


99 
. ZZ 


9A1 
ZOl . 


no 
. uu 


n 
u. 


1 A 
. 14 


1 
1 


8 1 

.Ol 


M099/I 

1V1UZZ4 


9 -1 
Z4. 


■ ol 


OU4. 


7^ 


n 
u 


O 9 
.UZ 


997 
ZZ 1 . 


'3 A 
. ol 


n 
u 


9 '3 
. Zo 


nqn 
ZoU. 


1 K 

. lo 


U- 


1 A 
. 14 


n 
u 


O 

.yo 


M0997 
1V1UZZ / 


OU. 


88 

on 


^I 7 
Ol i . 


. IO 


n 
-u 


on 
.uu 


Z44. 


7 '3 
. f O 


n 
u 


1 S 
. IO 


947 
Z4 ( . 


1 

. iy 


U. 


1 8 

. lo 


i 


"3("1 

. o u 


M09^G 

iviuzoy 


iy . 


Q O 
.OO 


907 

zy ( . 


Q '3 

. OO 


-U 


OA 

.uu 


997 
ZZ 1 . 


1 A 

. 10 


u 


1 '3 

. lo 


99A 
ZZD. 


07 

. y / 


U. 


1 9 

. 1 Z 


i 

i 


7 '3 

. / o 


iviuzyu 


99 


n i 


oiy. 


7 r ! 


n 
~u 


1 ("1 


99^ 
ZZO. 


'31 
Ol 


n 


1 9 
. 1Z 


99Q 

zzy. 


9-1 
. Z4 


n 
u. 


1 1 

. 1 1 


i 


AA 
.OO 


moqoo 

IVlUoUU 


ICS. 


ft 

■ OO 


977 
Z / / . 


qfi 
OU 


n 

-U 


! n 

. 1U 


90/1 
ZU4. 


1 H 
. lo 


n 
u 


1 1 
. 1 1 


907 
ZU I - 


1 O 
. 1U 


U- 


no 
uy 


i 


■3 7 
.O ( 


iviuozy 


9 "1 
Zl . 


'3 1 
.Ol 


ouy. 


"^0 

. oy 


-u 


• UO 


9*3^ 
ZOO. 


AO 

. oy 


u. 


1 8 
. IO 


9*37 
ZO i . 


97 
. Z ( 


u. 


1 1 

. 1 1 


o 
Z 


1 ^ 
• 10 


MOQBO 
IVlUooU 


"1 7 
1 / . 


08 

.Uo 


ouy. 


77 


("1 
-U 


OA 

.uu 


9/1 
Z4y. 


41 


u 


1 n 

. 1U 


9/1 C 
Z4o. 


Ol 


U. 


no 
. uy 


9 
Z 


O '3 

.Uo 


M0408 

1V1U4UO 


"1 7 


71 

. * i 


907 

zy t . 


90 

. zu 


n 


"1 '3 
. IO 


91 7 
Zl i . 


A '3 
.10 


n 
u 


OA 
.UO 


9 1 

ziy. 


AO 
. OU 


n 
u. 


UA 

. uo 


1 

1 


An 

.OU 


1V1U440 


9 '3 

ZO. 


.Uo 


o40. 


n*3 

. Uo 


-U 


OO 

.uy 


9 '3 % 
Zoo. 


■ OO 


n 
u 


on 
.zu 


9Q7 
Zo I - 


",A 
■ OO 


U- 


no 
,U9 


1 

1 


.0 1 


iviuo4y 


1 "1 


A 1 

1 


9"37 
ZO i . 


'3-1 
. 04 


n 

u, 


1 
.U 1 


1 OO 

iyu. 


'39 
. OZ 


n 
u 


1 

. iy 


1 80 

ioy. 


9 '3 
. ZO 


U. 


1 "3 

. IO 


1 

1 


07 
.U ( 


A /THAI A 
IVIUOIO 


1 '3 

lo. 


.Uo 


zoy. 


70 


-U 


o o 
.uy 


9 n - 3 

ZUo. 


. OZ 


u 


o o 
.uy 


onA 

ZUO. 


9-1 
. Z4 


U. 


1 1 

. 1 1 


1 


9 8 
.Zo 


1V1UOD4 


1U. 


.oy 


94A 
Z4D. 


"1 A 

IO 


n 
-u 


00 

.uy 


1 84 
104. 


1 8 

■lo 


n 


1 A 
. 14 


1 87 
lo 1 . 


1 r ! 

.10 


n 
u. 


. U i 


1 
1 


"3 A 
OO 


MO 791 
1V1U / Zl 


lo. 


■3 7 

.0 / 


97A 
Z /O. 


"1 9 
. 1Z 


n 

-U 


9 1 
. Z 4 


I - 3 

iyo. 


.40 


n 

-U 


o c: 
.Uo 


1 OQ 

iyo. 


AQ 

. Oo 


-U 


09 
.UZ 


9 
Z 


. lo 


M07A'i 
IVlu / DO 


1 '3 

io. 


A8 

■ Do 


9 '3 8 
ZOO. 


09 

. yz 


-u 


00 

.uy 


"I 70 

i / y. 


79 
. ( Z 


n 
u 


01 

.U4 


1 8*3 
loo. 


1 1 

.14 


U. 


n i 

. Ul 


1 

1 


70 


IVlUoOO 


"1 1 
14 


9 A 
.ZO 


9A4 
ZU4. 


88 

. OO 


n 
-u 


9 '3 
.ZO 


1 71 

1(1. 


00 


u. 


09 
.UZ 


1 7A 
1 / D. 


'30 

. oy 


n 

-u. 


1 
.Ul 


9 
Z 


0'3 

.uo 


Moons 


1 '3 

IO. 


A*} 
.10 


9A4 
ZU4. 


A 9 
. 1Z 


n 


9 1 
. Z 1 


1 7A 
I/O. 


(V3 

. uo 


n 
u 


09 
.UZ 


1 78 
I/O. 


1 s. 

. lo 


n 
u. 


n 1 

in 


1 

1 


A 8 
.UO 


M0G48 

iviuy4o 


Q 

y . 


9 '3 
. Z O 


9*37 
ZO ( 


98 
. ZO 


n 
~u 


09 
.UZ 


1 QH 

iyo. 


- r. 
. OO 


n 


1 1 
. 1 1 


1 07 

iy / . 


08 

. yo 


n 
u. 


1 
. 1U 


1 

1 


OA 

.yo 


MOQ^Q 

iviuyoy 


o 

O. 


1 "1 
.11 


999 

zzz. 


79 
. f Z 


n 
-u 


OA 

.uo 


1 A7 


71 
. i 1 


n 
u 


1 ^ 
■ 10 


1 71 
1/1. 


■1 7 
.4 / 


r\ 
u. 


("19 
. UZ 


9 
Z 


1 1 
. 1 1 


Mi m 7 

1V11U1 / 


O 
O. 


87 
.0 / 


9^4 
Z04. 


. OU 


-u 


. lo 


1 89 
loZ. 


^A 

. OU 


-u 


1 
.Ul 


18^ 
loO. 


A r i 
00 


u. 


n i 

. Ul 


1 

1 


. o y 


M1061 


7. 


.20 


206. 


25 


-0 


.15 


157. 


.89 





.02 


159. 


.55 


0. 


.04 


l 


.72 


M1071 


10 


.82 


240. 


.60 


-0 


.17 


157. 


.46 





.06 


159. 


.63 


0. 


.07 


l 


.26 


M1091 


10. 


.1(5 


243. 


14 


-0 


.21 


152. 


.81 





.01 


155. 


.36 


-0 


.02 


i 


.13 


M1167 


10. 


.21 


249. 


.04 


-0 


.17 


172. 


19 





.06 


173. 


.63 


0. 


.04 


i 


.56 


M1192 


6 


.05 


189. 


.01 


-0 


.07 


141. 


12 





.12 


144. 


.20 


0. 


.06 


l 


,16 


M1196 


10 


.74 


257. 


.84 


-0 


.19 


188. 


.89 


-0 


.01 


189. 


.65 


0. 


00 


2 


.02 


M1306 


9. 


.04 


245. 


04 


-0 


.21 


156. 


.97 





.01 


160. 


81 


0. 


.01 


1 


.52 


M1646 


7. 


.47 


224. 


.12 


-0 


.20 


161. 


.42 





.01 


164. 


.57 


-0 


.02 


1 


.80 


M1859 


6 


.27 


211. 


.94 


-0 


.13 


149. 


.00 





.08 


152. 


.64 


0. 


01 


1 


.41 


M2283 


4. 


.70 


187. 


.63 


-0 


.17 


126. 


.38 





.05 


128. 


99 


-0 


.05 


1 


.72 


M2665 


4. 


.40 


185. 


.07 


-0 


.22 


131. 


.70 


-0 


.03 


135. 


.47 


-0 


.02 


1 


.84 


M3431 


2 


.59 


175. 


.82 


-0 


.25 


119. 


.96 


-0 


.02 


125. 


09 


-0 


.10 


1 


.87 


M3852 


3. 


.6(5 


173. 


.24 


-0 


.23 


123. 


.50 


-0 


.09 


127. 


.22 


-0 


.15 


1 


.97 


M4323 


3 


.31 


169. 


72 


-0 


.27 


114. 


.39 


-0 


.09 


116. 


.01 


-0 


.08 


1 


.45 


M5014 


3 


.14 


163 


.71 


-0 


.30 


109. 


.17 


-0 


.12 


116. 


.20 


-0 


.15 


2 


.41 


M6782 


2 


.71 


158. 


.66 


-0 


.28 


112. 


11 


-0 


.10 


116. 


.33 


-0 


.1(5 


1 


.97 



Caon N., Capaccioli M., D'Onofrio M., 1993, MNRAS, 265, 
1013 

Capaccioli M., 1989, in H. G. Corwin Jr. & L. Bottinelli ed., 
World of Galaxies (Le Monde des Galaxies) Photometry 
of early- type galaxies and the R exp 1/4 law. pp 208-227 

Cappellari M., Bacon R., Bureau M., Damen M. C, Davies 
R. L., de Zeeuw P. T., Emsellem E., Falcon-Barroso J., 
Krajnovic D., Kuntschner H., McDermid R. M., Peletier 
R. F., Sarzi M., van den Bosch R. C. E., van de Ven G., 
2006, MNRAS, 366, 1126 

Cappellari M., Emsellem E., Bacon R., Bureau M., Davies 
R. L., de Zeeuw P. T., Falcon-Barroso J., Krajnovic D., 
Kuntschner H., McDermid R. M., Peletier R. F., Sarzi M., 



van den Bosch R. C. E., van de Ven G., 2007, MNRAS, 
379, 418 

Cassata P., Giavalisco M., Guo Y., Renzini A., Ferguson 
H., Koekemoer A. M., Salimbeni S., Scarlata C, Grogin 
N. A., Conselice C. J., Dahlen T., Lotz J. M., Dickinson 
M., Lin L., 2011, ApJ, 743, 96 

Churazov E., Forman W., Vikhlinin A., Tremaine S., Ger- 
hard O., Jones C, 2008, MNRAS, 388, 1062 

Churazov E., Tremaine S., Forman W., Gerhard O., Das 
P., Vikhlinin A., Jones C, Bohringer H., Gebhardt K., 
2010, MNRAS, 404, 1165 

Coccato L., Gerhard O., Arnaboldi M., 2012, ArXiv e- 
prints 



© 0000 RAS, MNRAS 000, 000-000 



16 Wu et al. 



Table 2. Fit parameters for stellar surface density profiles, ellipticities, and rotation parameters for the simulated galaxy sample: 
designations of the model galaxies (1st column), best fit parameters for cored-Sersic profiles: normalization parameter for central surface 
density, Jo, Sersic index n, core radius Ro, effectiv e radius R' e (2n d-5th columns); the 6th column is R e (i.e., the effective radius 
containing half the stellar mass within 10% r vir from lOser et alj|2010h : 7th-8th columns give ellipticities at 1 R c and 5 R c for the stellar 
mass distributions, 9th and 10th columns give rotation parameters A(5 R c ), and ratio v/a(5 R e )/?;/cr( R e ). 
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